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ABSTRACT 

We explore galaxy properties in general and properties of host galaxies of gamma- 
ray bursts (GRBs) in particular, using N-body/Eulerian hydrodynamic simulations 
and the stellar population synthesis model, Starburst99 , to infer observable properties. 
We identify simulated galaxies that have optical star formation rate (SFR) and SFR- 
to-luminosity ratio similar to those observed in a well-defined sample of ten host 
galaxies. Each of the numerical counterparts are found in catalogs at the same redshifts 
as the observed hosts. The counterparts are found to be low-mass galaxies, with low 
mass-to-light ratio, recent epoch of formation, and high ratio between the SFR and the 
average of the SFR. When compared to the overall galaxy population, they have colors 
much bluer than the high-mass star-forming galaxy population. Although their SFRs 
span a range of values, the specific rates of the numerical counterparts are equal to or 
higher than the median values estimated at the different redshifts. We also emphasize 
the strong relationships between the specific star formation rate (SFR) and quantities 
known to reflect the star formation history of galaxies, i.e. color and mass-to-light 
ratio: At intermediate redshift, the faintest and bluest galaxies are also the objects 
with the highest specific rates. These results suggest that GRB host galaxies are likely 
to be drawn from the high specific SFR sub-population of galaxies, rather than the 
high SFR galaxy population. Finally, as indicated by our catalogs, in an extended 
sample, the majority of GRB host galaxies is expected to have specific SFRs higher 
than found in the magnitude-limited sample studied here. 

Key words: cosmology: large-scale structure of Universe - galaxies: formation - 
galaxies: evolution - gamma rays: bursts 



1 INTRODUCTION 

Galaxies with a wide variety of properties are observed in 
the universe and galaxy sub-populations contribute in a 
different way to the overall galaxy properties at different 
redshifts. Important questions in present day cosmology in- 
clude investigations into how the different sub-populations 
can be characterized, how sub-populations formed at high 
redshift evolve into galaxies in the local universe and how 
these populations contribute to various properties of the 
overall galaxy population. In this paper we use a numer- 
ical ap proach to extend o ur investigation of galaxy prop- 
erties (|Courtv et al.l |2004| . hereafter referred to as Paper 
I) that focused on the specific star formation rate. While 
the star formation rate is a 'snapshot' of the stellar activ- 
ity in a galaxy, the specific rate is an indicator of how the 
galaxy forms its stellar mass relative to the total mass that 
has been assembled through its entire lifetime, via mergers 



and/or transformation of new accreted gas. The specific rate 
should therefore give some insight into the star formation 
history of galaxies, as do other properties, such as colors. 
The specific SFR has been estimated in a number of ob- 
servational studies: iGuzman et all l|l997l) compare the spe- 
cific SF Rs of compact blue galaxie s with other galaxy popu- 
lations, iBrinchmann et al. (12004 ) discus s the star- forming 
galaxies of the SPSS. iBell et alj (|2005h combine infrared 
dat a from the Sp itzcr Space Telescope with optical data, 
and iFeulner et alT i2005~ estimate the specific SFR up to 
2 = 5. 

As in Paper I, we here focus on the properties of a 
particular population of galaxies, namely the host galax- 
ies of long-duration gamma-ray bursts. The nature of the 
hosts and their evolution with redshift are still open ques- 
tions. The emphasis over the last few years has been on 
the host galaxies of long-duration GRBs, which are seen 
as a powerful tracer of massive star formation in the uni- 



2 Courty, Bjdrnsson & Gudmundsson 



verse. It is now well established that at least some long- 
duration GRBs occur at the death o f massive stars (e.g. 
iHiorth et all l2003i : IStanek et all 120031 ) and in a cosmolog- 
ical context massive stars are very short-lived. Because of 
their extreme brightness, long-duration GRBs are an effec- 
tive way of locating distant galaxies, most of which are so 
faint that they would go undetected in galaxy surveys. In 
fact, host gala xies fainter than ma gnitude 29 have already 
been detected (|jaunsen et al.l l2003). In addition, GRBs have 
been detected out to a redshift higher than z — 6, and will 
likely be detected to even higher redshifts. Rather detailed 
studies of the local star forming r e gions in the hosts is there- 
fore possible (|Berger et all l|2005l ). lBerger et all l|2005l) ). The 
most distant burst to date, GRB 0509 04 at z = 6.29 
(|Tagliaferri et al.ll2005l ; iKawai et "all 120051 ). has already re- 
vealed a number of i nteresting pr operties of the interstellar 
medium of the host (Totani 2005), while information about 
galaxy formation and evolution is yet to be explored in de- 
tail, mostly awaiting the increase of the currently modest 
sample size (around 50 hosts). 

;,From individual studies of host galaxies of GRBs 
(|Fruchter et al J 19991 : iFvnbo et ai1l2003l ) , comparison of host 
samples with other sou rces detected in various deep surveys 
|Le Floc'h et al.l l2003h . and statistical stellar population 
synthesis of the opt i cal and near-infra-red host properties 
JSokolov et all l200ll ; ICharv et all 120021 : IChristensen et all 
120041 '). indications are that host galaxies have particular 
characteristics: These galaxies tend to be optically sub- 
luminous, low-mass, blue, star-bursting galaxies, with young 
stellar populations, a modest activity of optical star forma- 
tion, and perhaps low-metalicity and modest amount of dust 
obscuration, although this la s t feat ure still needs to be firmly 
established. iLe Floc'h etaH (|2003l ) compare a large sample 
of host galaxies of GRBs, observed in the near infra-red, 
with various galaxy surveys and find that the observed K 
and R— band magnitudes of the hosts are comparable to the 
field sources selected in optical/N-IR deep surveys, but dif- 
fer significantly from luminous and dusty starburst galaxies 
observed with ISO and SCUBA. Moreover they show blue 
R — K colors typical of the faint blue galaxy population in 
the field at z = 1. Also, SCUBA su b-millimeter observa tions 
of GRB host galaxies performed bv lSmith et al.l (|2005l ) sug- 
gest that most hosts are not luminous dusty s tar-forming 
galaxies. Furthermore, IChristensen et all l|2004h have stud- 
ied a magnitude-limited sample of hosts and estimate the 
ratio between the rest-frame UV star formation rate and 
the host optical luminosity. They suggest that the hosts are 
similar to those HDF galaxies that have the highest SFR- 
to-luminosity ratios. 

In Paper I, we showed that among the population still 
actively forming stars at low redshift, the high-mass galaxies 
have much lower specific SFR than the low-mass galaxies. 
The non star-forming galaxies, that span the whole galaxy 
mass range, are old galaxies while most of the stellar pop- 
ulations in the high specific SFR galaxies formed recently. 
At high redshift the trend of increasing specific SFR with 
decreasing galaxy mass is also seen, but an interesting point 
is that the cosmological evolution is much stronger for the 
high-mass than the low-mass galaxies. These trends agree 
in general well with the aforementioned observational esti- 
mates, to the extent that we concentrate on the qualitative 
behavior of the specific SFR. Although not based on ob- 



servable properties, the results of Paper I for the specific 
SFR suggest that a sub-population of faint and blue galax- 
ies, some of the characteristics of the GRB hosts, are likely 
to belong to the high specific SFR galaxy population, rather 
than the high-SFR population. 

In this paper, we extend the discussion of the properties 
of GRB host galaxies in Paper I, by combining the results of 
the same simulations as in that paper with t he stellar pop- 
ulatio n synthesis (SPS) code, Starburst99 (|Vazquez et al.l 
2005), to infer observable properties. SPS codes allow us to 
compute the spectral energy distributions of the simulated 
galaxies at different redshifts. We identify simulated galaxies 
that have both similar rest-frame ultraviolet SFRs and ra- 
tios between this SFR and the B— band lu minosity as the ten 
observ ed hosts of a well-defined sample (|Christensen et all 
(|2004l ). hereafter Chr04). The numerical counterparts of 
these observed hosts are characterized by a variety of prop- 
erties, estimated either directly from the simulation or from 
the computation of the SEDs. In particular, for each galaxy 
we determine the mass, the ratio between the SFR and 
the average of the SFR, SFR* /(SFR), the specific SFR, 
the epoch of formation, the R — K color, and the mass- 
to-light ratio. The properties of the counterparts are then 
compared to the overall galaxy population, that is charac- 
terized through the close relationships between the specific 
SFR and the color index and mass-to-light ratio. This com- 
parison is, however, limited by the fact that the observed 
sample still only includes 10 hosts, spanning a large redshift 
range. 

Although fairly small, the sample of Chr04 is the only 
available homogeneous sample that estimates the ratio be- 
tween the optical SFR and the luminosity. Information on 
other hosts does exist i n the literature and other studies, e.g. 
ISollerman et~"ai1 (|2005l ). use different star formation rate es- 
timators such as SFRn a and SFRo u , rather than the UV- 
based indicator adopted in Chr04, but only for a couple of 
hosts. A comprehensive study of a large host sample using 
all available information awaits future studies. In Section 
[3] we also briefly discuss the counterparts of GRB 000911 
a nd 030329, whose SF Ru y and Mb values are available 
in lMasetti et all (120051 ) and lGorosabel et al.1 (|2005l ). respec- 
tively. 

The paper is organized as follows: In Section [5] we de- 
scribe briefly the simulation and use of the SPS code. Sec- 
tion [3] contains the main results of the paper. It starts by 
discussing the observed host sample we use (section I3.1|l 
and the procedure used to identify the numerical counter- 
parts (section T3.2|l . We discuss the observational properties 
of the counterparts in section 13.31 and compare their prop- 
erties with those of the galaxy population in section 13.41 
Section [4] concludes the paper. 



2 NUMERICAL PROCEDURE 

We repeat the simulation used in Paper I in order to 
obtain galaxy catalogs at the same redshifts as that of 
the observed GRB hosts in the sample of Chr04. We 
briefly recall the numerical method and we refer to Pa- 
per I for details regarding the simulation. The three dimen- 
sional N-body/hydrodynamical code couples a PM scheme 
for computing gravitational forces with an Eulerian ap- 
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proach for solving the hydrodynamical equations. The dom- 
inant processes relevant for galaxy formation are included: 
Gravitation, hydrodynamical shocks and radiative cool- 
ing processes. Collisional ionization equilibrium is not as- 
sumed and the cooling rates are explicitly computed from 
the evolution of a primordial composition hydrogen-helium 
plasma. The cosmological scenario adopted is a A— cold 
dark matter model with the following parameters: Ho = 
70 km s" 1 Mpc -1 , Q K = 0, Q. m = 0.3, Q. A = 0.7, Q. b = 
0.02/i" 2 with h = Ho/100, cr 8 = 0.91. The comoving size 
of the computational volume is 32 /i~ 1 Mpc and the simu- 
lation has 256 3 dark matter particles and an equal number 
of grid cells. Galaxy formation is introduced using a phe- 
nomenological approach. At each time step, in cells whose 
gas satisfies given criteria, a fraction of the gas is turned into 
a stellar particle of mass m, and epoch of formation a,, the 
scale-factor corresponding to the cosmic time, t*. The crite- 
ria are the following: The cooling time must be less than the 
dynamical time, i coo i < td yn ', the baryonic density contrast, 
(1 + 8b), must be higher than a threshold (1 + 5b)s = 5.5; 
the gas must be in a converging flow, V ■ v < 0; and the 
size of the cell must be less than the Jean's length, given by 
Aj = Csin/Gp) 1 ^ 2 . The mass m* is given by mB(to)At/t,, 
where At is the timestep, the characteristic time is taken 
to be i, = max(£dyn, 10 8 yr), and tub (to) is the baryonic 
mass enclosed within the grid cell. Galaxy-like objects are 
then defined, at any redshift, by grouping the stellar parti- 
cles with a friend-of-friend algorithm. A simulated galaxy is 
therefore a collection of stellar particles of different masses 
formed at different epochs. 

At any redshift, galaxies are characterized by their 
mass, M, defined as the sum of the mass of all the stel- 
lar particles the galaxy is composed of; an epoch of forma- 
tion, defined to be the mass-weighted average of the epoch 
of formation of all its stellar particles; and a star forma- 
tion rate, SFR*, defined as the amount of stellar material 
formed in the previous 10 8 yr. We also estimate the specific 
star formation rate, e = log(10 n yr SFR* /M), a quantity 
that measures the efficiency of the conversion of the gas into 
stellar material relative to the galaxy mass. We have, as in 
Paper I, only considered galaxies with a mass higher than 
5 ■ 10 s Mq in the simulated catalogs. We refer to Paper I for 
an illustration of the different sub-populations of galaxies 
obtained. 

Considering each stellar particle as a homogeneous stel- 
lar population, the total spectral energy distribution (SED) 
of a galaxy can be computed from synthesis codes. Evolu- 
tionary synthesis codes combine stellar evolutionary theory 
to describe the time evolution of model stars, stellar atmo- 
sphere theory to transform quantities from the theoretical 
space to the observational one, and the stellar birth rate, 
giving the number of stars with a given initial mass formed 
at a given time. The stellar birth rate comes from the star 
formation history, which gives the number of stars born in 
a given time, and the initial mass function (IMF), which 
give s the relative number o f star s born as a function of mass 
(see[c crvmo fc Luridianal (|2005h for a discussion regarding 
the uncertainties of the synthesis codes). In our procedure, 
the star formation history of a galaxy comes directly from 
the formation epoch and mass of all its stellar populations 
as recorded in the simulation. We use the stellar popu la- 
tion synthesis model, Starburst99 jLeitherer et al.lll999l ) in 



its latest version l| Vazquez et al.ll2005h . to derive observable 
properties. The main change in this version is the introduc- 
tion of the Padova stellar evolutionary tracks, allowing the 
computation of stellar evolution for old and low-mass stars 
as well as high-mass stars. 

We start by defining a SED template: we consider a sin- 
gle stellar population of 10 6 Mq, assumed to form instan- 
taneously and evolve passively over a maximal time of 15 
Gyr. The SED is computed at 1221 points between 91 A and 
160 [im. We select in Starburst99 the "Padova AGB" evo- 
lutionary tracks, selection of the 1992 — 1994 Padova tracks 
with thermally pulsing AGB stars added, with a metalic- 
ity of Z = 0.004 = 0.2Z e , where Z e ~ 0.02. We choose 
a Salpeter IMF, dP/dm = m -2,35 , with low and high-mass 
cut-offs at 0.1 and 100 Mq, respectively. The other input pa- 
rameters adopted (supernova cut-off mass, black hole cut-off 
mass, wind model, interpolation in mass method, model at- 
mosphere), are the standard ones. We refer the reader to the 
Starburst99 websitcQ for a full explanation. 

The total SED of a galaxy at a given redshift is the 
sum of the SEDs of all the stellar populations present in the 
galaxy. The SED of a stellar particle, at a given redshift, is 
the SED template evolved on the time t(z) — t*, in order to 
take cosmological expansion into account, and scaled to the 
mass of the stellar population with the factor m» / '(10 6 'Mq). 
To properly include the contribution of those stellar particles 
that formed at the simulation output time, t(z), the SED is 
evolved over 10 4 yr. Note that this procedure is applied to 
the simulated galaxies after the simulation run, and not dur- 
ing it. Nevertheless, the inferred SED of a galaxy accounts 
for the formation all along the simulation of different stellar 
populations. A caveat of the procedure is that the simulation 
does not take into account any evolution of the gas metal- 
icity, and our derived observable properties are based on a 
single template. Due to the higher UV-luminosities of low- 
met alicity stars, assuming higher metalicity tracks would de- 
crease the total UV-luminosity. However, in the catalog at 
z = 1 almost all of the non star-forming galaxies have lumi- 
nosities more than twice the luminosity computed with so- 
lar metalicity tracks, whereas there are only 19 star-forming 
galaxies with such a large luminosity. 

Star formation rates can be derived from the UV contin- 
uum or the Ha luminosity using appropriate proportionality 
coefficients. These coefficients are derived from stellar pop- 
ulation synthesis models, involving an initial mass function 
(IMF) and star formation history, by looking for a propor- 
tionality relation assuming that the star formation rate is 
constant in the later phases on given timescales. These are 
roughly 10 7 yr for Ha and 10 8 yr for the UV. These coef- 
ficients therefore depend on the IMF and its s lope, the up- 
per a nd lower mass cut-off, and the metalicity (jBoselli et alj 
l200lT ). Since one of the purposes of this paper is to find the 
numerical counterparts of observed host galaxies, we use the 
same indicator of the star formation rate as Chr04, and es- 
timate it from the UV continuum at 2800 A, although the 
SFRs of the two hosts at the highest redshifts are estimated 
using the UV continuum at 1500 A in Chr04. The calibra- 
tion factor used to convert the luminosity at 2800 A into 
SFR is 1.4 x 10" 28 , in units of (Mq yr _1 )(erg s _1 Hz -1 ) -1 
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SFR* (M e yr l ) 

Figure 1. Comparison of the star formation rate SFR* with 
SFRuv obtained from the synthetic spectrum for galaxies at 
z = 1. Non-star forming galaxies are for display purposes plot- 
ted at SFR* = 0.02 M Q yr -1 . The diagonal dashed line shows 
SFR UV = SFR*. 



l|Kennicuttlll99ct ). In our procedure the luminosity at 2800 
A is obtained by averaging the spectrum over 20 A around 
2800 A. Note that, although our spectral energy distribu- 
tions are determined at sub-solar metalicity, we adopt the 
usual calibration to convert the UV luminosity into SFR 
that is estimated at solar metalicity. At lower metalicity the 
calibration increases due to the higher luminosity of low- 
metalicity stars. The SFR of low-metalicity galaxies when 
estimated with a calibration at solar metalicity is thus over- 
estimated. The evolution of the calibr ation associated with 
diff erent SFR indicators is estimated in Sullivan et alj (|200ll ) 
and lBicker fc Fritze-v. Alvenslebe for instance. At 

Z — 0.2Zq, the latter paper found that the SFR is overes- 
timated by a factor 1.4, compared to our calibration. 

In Figure [T] we compare the SFR* , computed directly 
from the simulation, with the SFR estimated from the syn- 
thetic spectrum of the simulated galaxies at z = 1. Non- 
star forming galaxies (i.e. with SFR*—0) are for display 
purposes plotted at SFR* = 0.02 M yr -1 . At this red- 
shift, the galaxy population includes 1927 objects, with 1148 
star-forming galaxies (as expected, the number of non star- 
forming galaxies is larger at z = 0). The diagonal dashed 
line shows SFRuv = SFR* . Despite the scatter in the dia- 
gram, it is clear that the two different estimates are of the 
same order. The estimate from the simulations is a factor of 
about 2 lower than that obtained from the synthetic spectra. 

Figure [2] shows the B — and .ft'— band magnitudes as 
functions of the galaxy mass for the simulated catalog at 
2 = 1. The B and K— band magnitudes are computed using 
the "Buser's S3" filter and the "IR K filter + Palomar 200 
IR detectors + atmosphere 0.57" from the Galaxev pack- 
age (|Bruzual fc Cha riot 2003). Different symbols distinguish 
star-forming and non star-forming galaxies. Note that both 
are strongly correlated with the mass. There is much less dis- 
persion in Mk since the B— band luminosities are dominated 
by the flux from young massive stars whereas Lk is a better 
tracer of the overall stellar population. The corresponding 
plots at z — would show the same trend although quantita- 
tively different. For instance, the Mk — mass relation shows 
slightly more dispersion and galaxies at M K = —24.2, the 




1 


10 9 


10 10 


10 11 


10 12 


■* 

CM 
1 




Mass (M Q ) 














CvJ 

cv 

1 










(mag) 
-20 










m 

S CD 
1 










CO 

' — 1 

1 








z=l 




10 9 


10 10 


10 11 


10 12 



Mass (M s ) 

Figure 2. The mass-magnitude diagrams (upper panel: K— band, 
lower-panel: B— band) for the simulated galaxies at z = 1. In both 
panels, dots denote star-forming galaxies while the small dots 
denote the non star-forming galaxies. The estim ate of M*^ = 
—24. 7 around 2 = 1 from the K20 galaxy survey l lPozzetti et al.1 
l2003t ) is marked by the dotted horizontal line. 



ch aracteristic magn itude of the galaxy luminosity function 
in lCole et all (|200ll 'l. have a mass of ~ 10 11 M . 

In the following sections we consider other quantities to 
characterize galaxies. The lifetime of a galaxy is estimated 
at a given redshift, t(z) —tf orm , where tf orm is the epoch of 
the formation of the first stellar population. From this the 
average of the star formation rate, {SFR), at a given red- 
shift is determined as the amount of stellar material formed 
over the lifetime of the galaxy. The SFR* is compared with 
the average through the ratio SFR* /{SFR}, that may be 
seen as analogous to t he so-called birth-rate parameter, b, 
l|Kennicutt et al.lll994l ). The mass-to-light ratios in the B 
and A'— bands are estimated using the solar magnitudes of 
5.45 and 3.3, respectively, and the color R — K using the 
"Cousins R" filter from the Galaxev package. Finally, the 
ratio between the rest-frame ultraviolet SFR and the lumi- 
nosity, SFRuv /{Lb / L* B ) , uses the characteristic magnitude 
or magnitude at the break of the local galaxy luminosity 
function, M B = —21. This ratio is sometimes called the 
specific SFR as e.g., in Chr04, but we will refer to it by the 
term SFR-to-luminosity ratio and reserve the term specific 
SFR for the ratio between the SFR and the galaxy mass, as 
is customary in the literature. 
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Table 1. Properties of the numerical counterparts corresponding to the GRB hosts listed in column one. Following the redshift of the 
simulated catalog in col. 2, the percentage errors in SFRuv and the SFR-to-luminosity ratio, AX and AY defined in the text, are given 
in cols. 3 and 4. Column 5 shows the mass of the counterpart, while col. 6 shows the epoch of its first formed stellar population, tf orm 
(see Section[2]for definition) relative to the Hubble time. The mass-to-light ratios in the K and B— bands, respectively, are given in cols. 
7 and 8, followed by the corresponding absolute magnitudes. Then comes the color index, R — K, and the next to last column lists the 
R(AB) apparent magnitude. Finally, a reference number for each host is given in the last column. 



GRB 


z 


AX 


Ay 


M (M ) 


1 — tforrn/^H(z) 


M/L K 


M/L B 


M K 


M B 


R-K 


R(AB) 


# 


000926 


2.036 


0.7 


0.5 


1.16 x 10 10 


0.76 


0.36 


0.36 


-22.97 


-20.83 


2.61 


24.41 


1 


990123 


1.592 


5.6 


3.2 


2.41 x 10 10 


0.81 


0.47 


0.59 


-23.47 


-21.08 


3.50 


24.27 


2 


000418 


1.113 


5.6 


8.8 


8.27 x 10 9 


0.86 


0.35 


0.32 


-22.63 


-20.57 


2.47 


23.24 


3 


980703 


0.963 


1.3 


0.6 


3.68 x 10 10 


0.88 


0.52 


0.66 


-23.82 


-21.41 


2.87 


22.18 


4 


000210 


0.843 


6.5 


5.7 


5.16 x 10 9 


0.89 


0.47 


0.54 


-21.79 


-19.50 


2.61 


23.64 


5 


970508 


0.832 


3.7 


0.3 


7.20 x 10 8 


0.89 


0.30 


0.25 


-20.14 


-18.18 


2.22 


24.84 


6 


991208 


0.704 


1.2 


0.9 


1.64 x 10 9 


0.90 


0.41 


0.45 


-20.71 


-18.46 


2.54 


24.24 


7 


970228 


0.695 


10.8 


0.5 


1.60 x 10 9 


0.90 


0.49 


0.56 


-20.49 


-18.19 


2.51 


24.41 


8 


010921 


0.449 


1.9 


1.1 


8.97 x 10 9 


0.92 


0.68 


1.01 


-22.01 


-19.42 


2.58 


21.93 


9 


990712 


0.429 


3.3 


2.1 


1.42 x 10 9 


0.92 


0.38 


0.36 


-20.62 


-18.55 


2.07 


22.71 
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3 THE NUMERICAL COUNTERPARTS OF 
OBSERVED GRB HOST GALAXIES 

In this section we attempt to identify numerical counter- 
parts to the host galaxies of the Ch04 sample, explore their 
properties, and compare the host candidate population with 
the overall galaxy population. We first summarize the rel- 
evant properties of the observed sample and our approach 
to identifying the numerical counterparts. We then discuss 
their observed properties. 



3.1 The observed sample 

The sample presented in Chr04 is homogeneous and con- 
sists of 10 GRB hosts with redshifts between z = 2.037 and 
2 = 0.433. The sample is magnitude limited (R < 25.3) and 
ensures that hosts are bright enough to make multi-color 
photometry possible. The first column in Table [1] lists the 
hosts under consideration. The observationally determined 
values of SFRuv and SFRuv I (Lb / ' L* B ) are given in Table 
4 and 5 of Chr04 and we note that the SFRuv entering the 
SFR-to-luminosity ratio is not corrected for internal extinc- 
tion, generally found to be moderate or low. The SFR spans 
a wide range (between 0.8 and 13 Mq yr _1 ), but hosts with 
redshifts higher than 0.9 have the highest SFRuv, more 
than ~ 6 Mq yr _1 , whereas hosts at lower redshifts have 
SFR UV < 2.5 M yr _1 . 

Various observational studies of individual hosts of the 
Chr04 sample re sult in slightly d i fferen t estimates of the 
same properties (IGorosabel et al.l (120031) for GRB 000418 , 
IGorosabel etafl (|2003h for GRB 000210. iBloom et all (|l999T ) 
for GRB 990123, for instance). We have not attempted 
to collect all the data existing in the literature, and have 
only considered the SFRuv and SFRuv / {Lb / L* B ) given 
in Chr04. All t he ho sts in Chr04 are also discussed in 
iLe Floc'h etafl (|2003h . giving either their observed R and 
absolute B— band magnitudes or their observed K magni- 
tudes and colors, R — K. In addition to the Chr04 sam- 
ple, we also consider two hosts for which the SFRs based 
on the rest-frame ultraviolet flux as well as the B— band 
magnitudes, are available. We can ex amine such hosts in a 
similar way as for the Chr04 sample. iMasetti et al.l (|2005l ) 
discuss the properties of the host galaxy of GRB 000911 at 



2 = 1.06 and derive a B— band magnitude of around —18.4, 
giving a Lb/L* b = 0.09 with our M B (see section [5}, and 
an extinction-corrected SFRuv of 2.7 Mq yr -1 . Assuming 
an extincted SFRuv to be a factor of 2 lowei0, gives a 
SFR-to-luminosity ratio around 15. The second h ost is that 
of GRB 030329 at z — 0.168. IGorosabel et al" | l|2005l ) de- 
rive a SFRuv of 0.17 M yr" 1 and M B = -16.5, giving 
a SFRuv / {Lb / L* B ) = 10.6. Note that this value is compa- 
rable to the SFR-to-luminosity ratios derived in the Chr04 
sample. The counterparts of these two hosts are discussed 
near the end of this section. 



3.2 Counterpart Identification 

We generated catalogs of galaxies at the same redshifts as 
the observed hosts, and looked in each catalog for simulated 
galaxies that have both SFRuv and SFRuv /(Lb/ L* B ) 
nearest to the corresponding values for each observed host. 
We are able to find a numerical counterpart to each of 
the 10 observed host of the Chr04 sample, although de- 
partures from the observationally inferred values may be 
large in some cases. The percentage errors in Table Q] quan- 
tify how close from these values the counterpart is found. 
They are defined by AX = \SFR°u b v - SFRuy 1 \/SFR u b v 
and AY = \ef s - el um \/e° L bs . In this last definition, e L de- 
notes the SFR-to-luminosity ratio. The largest departures, 
when at least one of these quantities is larger than 5%, are 
seen for counterparts #3, 5 and 8. The first two are noted 
in Chr04 as having the least acceptable fits to their spec- 
tral energy distributions The host of GRB 0004 18 (#3) is 
extensively discussed in IGorosabel et all (|2003T) who ana- 
lyze its spectral energy distribution using a variety of syn- 
thetic spectral templates. The adop ted SFR does of course 
depend on the adopted SED. Also, IGorosabel et all (I2003T ) 
find the host galaxy of GRB 000210 (#5) to be brighter 
than the Chr04 estimate: They find that the B— band mag- 



2 The host of GRB 991208 has roughly similar B-band mag- 
nitude and Ay and has a ratio of un-extincted/extincted SFR 
within a factor of 2 (Chr04). 
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Figure 3. SFR-to-luminosity ratio versus SFRuv (upper panel) 
and specific SFR versus SFR* (lower panel) for the star-forming 
galaxy population at z = 0.832 (dots). The dotted curves in the 
upper panel denote constant Ms-values from —18 to —23 (left 
to right). In the lower panel the dotted lines indicate constant 
mass of 10 9 , 10 10 and 10 11 M© (left to right). Although at dif- 
ferent redshifts, the numerical counterparts of the observed host 
galaxies (listed in Table [TJ are displayed in both panels (solid 
symbols): In decreasing order of SFR UV , GRB 980703, 000926, 
000418, 990123, 000210, 010921, 990712, 970508, 991208, 970228 
(the order is the same in the bottom panel with SFR*. except 
that the host of GRB 970228 is to the left of GRB 991208). In the 
upper panel, for each counterpart a short line segment indicates 
where the corresponding observed host galaxy would be found. In 
some cases the difference is too small for the line segment to be 
visible. The two diamonds denote the counterparts at the highest 
redshifts. 



rate SFR*, the specific SFR, the galaxy mass, the ratio 
SFR* /{SFR), the formation epoch of the first stellar pop- 
ulation to form in the galaxy, and the epoch of formation 
of the galaxy as defined in section [2] From the computation 
of the SEDs we determine the star formation rate SFRuv, 
the SFR-to-luminosity ratio, the magnitudes in the B and 
K— bands, the apparent R magnitude and the color, R — K. 
The mass-to-light ratios in the B and .ft'— bands are then 
a combination of the primary simulation outputs and the 
observable properties. 

3.3 'Observed' properties of the numerical 
counterparts 

From Table [1] we note the following: Low-redshift hosts 
(z < 0.9) are low B— band luminosity galaxies, Mb > —20. 
We also see that the counterparts are fainter than —22 in 
the K— band at z < 0.9. The 10 counterparts are low-mass 
galaxies with M < 4 x 10 10 Mq (the highest mass in the 
simulation at z = 1, being 7.4 x 10 11 M Q ). The B-band 
magnitudes and mass of the counterparts can also be seen 
in Fig. 13 discussed below, where in the upper panel the dot- 
ted curves denote values of constant Mb = —18 to —23 
(from left to right) and in the lower panel the dashed lines 
indicate values of constant mass of 10 to 10 1 Mq (from left 
to right). The counterparts have Lb/L* b between 0.074 and 
1.46 and L K /L* K between 0.028 and 0.84 (with M* K = -24), 
making most of them sub-luminous galaxies. The mass-to- 
light ratios in Table [T] are relatively small and more typical 
of late- type and dwarf galaxies than large spirals or elliptical 
galaxies. 

The counterparts have apparent i?-band magnitudes 
between 22 and 24.4 and comparison bet ween those and 
the N -IR observational data reported in iLe Floc'h et al.l 
(2003) for the same hosts (here referred to as the common 
sample) shows that the candidate hosts #6 and #9 are 
the faintest and the brightest, respectively, both of our 
sample and of the common sample. The apparent i?-band 
magnitudes of the cou nterparts are, howeve r, generally 
slightly brighter than in lLe Floc'h et all l|2003l ). The bluest 
of the 9 counterparts (#10) for which the R — K colors of 
the co rresponding hosts were estimated in ILe Floc'h et al.l 
(2003), is also the bluest of these 9 observed hosts. The 
observed host with the highest color index corresponds 
the counterpart with the second highest color index 
(after the least blue #2). Globally, we do find that 9 coun- 
terparts have R — K colors between 2 and 2.9 and 3.5 for #2. 



nitude is —20.16, giving, with their SFRuv =2.1 M Q yr , 
a SFR uv /(L b /L*b) = 4.5 (10.7 in ChrO4).0 

Each numerical counterpart is characterized by a num- 
ber of properties that are either tabulated in Table [1] or dis- 
played in Fig. [3] and [5] The properties that are estimated 
directly from the simulation output are the star formation 



3 Instead of adopting the nearest counterpart in each catalog, we 
could have looked for the second best closest counterparts or the 
counterparts within a given error box around each observed host. 
Because of the uncertainties in the observationally determined 
values and the small number of observed hosts in the sample we 
consider here, we do think that searching only for the closest 
counterpart is the most reasonable thing to do. 



In Fig. [3] we superpose the numerical counterparts on a 
plot showing the simulated star-forming galaxy population 
at z = 0.832 in a diagram of SFR-to-luminosity ratio versus 
SFRuv (upper panel) and in a diagram of specific SFR ver- 
sus SFR* (bottom panel). We have chosen z = 0.832 as it is 
roughly the median redshift of the sample. Rigorously speak- 
ing, only the counterpart of the host at z = 0.832 should be 
compared directly with the simulated catalog, since galaxy 
properties evolve with redshift. It may nevertheless be useful 
to plot the whole sample on a single diagram. We have used 
different symbols (diamonds) to mark the highest two red- 
shifts, and will return to the issue of evolution later in this 
section. Each counterpart in the top panel is accompanied 
by a straight line, sometimes smaller than the symbol, that 
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indicates how well the properties of that numerical coun- 
terpart match those of the observed host (the difficulty in 
obtaining a good counterpart for host #3 is clearly seen). As 
noted above, the hosts have various SFRs, be it either the 
SFRuv or the SFR* , but they have high SFR-to- luminosity 
ratio and high specific rate, as they all lie in the upper half 
of both panels in Fig. [3] The specific SFR and the SFR-to- 
luminosity ratio would be even better correlated if the latter 
quantity is based on the A"— band luminosity, according to 
the magnitude-mass diagrams in Fig. [2] 

3.4 Comparison with the overall galaxy 
population 

Even if some aspects of the simulation may be incomplete, it 
provides us with catalogs of galaxies at various redshifts in 
which different galaxy populations may be distinguished. Af- 
ter defining a particular sub-population of galaxies e.g. sim- 
ilar to the host galaxies of GRBs, an important issue is 
to compare this particular population to the overall galaxy 
population. Such a comparison is in the present case limited 
by the fact that the observed sample includes only 10 hosts, 
spanning a wide redshift range where cosmological evolution 
cannot be neglected. 



3.4-1 The relationship between the specific SFR and the 
color index 

We first conduct a qualitative comparison between the coun- 
terparts and the overall galaxy population at z = 1, focus- 
ing on the relationship between the specific SFR and the 
mass-to-light ratio and color index. By presenting observ- 
able properties at this redshift we follow-up on and extend 
the results and discussion of Paper I. 

Figure [4] presents the mass-to- light ratios and color in- 
dex of the whole galaxy population distinguished accord- 
ing to the specific rate (e < 1, 1 < e < 1.3 and e > 1.3). 
We also plot the low-mass (M < 5 ■ 10 10 M ) and high- 
mass (M > 5 • 10 10 M©), non-star forming galaxies. We use 
the same threshold, e = 1.3, as in Paper I. It corresponds 
roughly to the peak of the probability density function of 
the specific SFR at z — 1. In addition, galaxies with specific 
rates below or above this value contribute about equally to 
the total star formation rate at that redshift. The top panel 
in Fig. U shows the large variation of the mass-to- light ratio 
with respect to the galaxy mass, but there is a clear cor- 
relation with the specific SFR (i.e. e). High- mass non star- 
forming galaxies have large M/Lk, up to about 1, whereas 
the minimal value for the low-mass high-specific rate galax- 
ies is ~ 0.15. The mass-to-light ratio in the B— band, M/ Lb, 
ranges between ~0.09 and ~3, for the overall population. 
Recalling the results in Table [1] we note that the major- 
ity of the counterparts have \og(M / Lk) between —0.5 and 
—0.3, typical of our low-mass, high specific SFR galaxies, 
although it should be kept in mind that Fig. [4] is plotted at 
a single redshift (z = 1). 

The color-magnitude diagram in the bottom panel in 
Fig. [4] shows a large variation of the color index R — K, 
from 1.5 up to 4.7, consistent with that seen in the recent 
K20 galaxy survey (|Pozzetti et alj [20031 ). The high-mass, 
non star-forming galaxies, that are also old objects (see Pa- 
per I) have the highest color index, close to the typical value 



o 
s 



■■ ■ ■■ ■" " ' 



icr 



10 1 



Mass (M„) 



10 1 



10 1 



I 

K 



.... 



. * t *. ■ 



-IE 



-20 



-22 



-24 



-26 



M K (mag) 



Figure 4. Mass-to-light ratio in solar units in the K— band as 
a function of mass (upper panel) and color-magnitude diagram 
(bottom panel) for the simulated galaxy population at z = 1. 
For clarity, only 1000 objects randomly selected from the catalog, 
are plotted. The dots indicate non star-forming, low-mass (M < 
5 • 10 10 Mq) galaxies; the filled circles non star-forming, high- 
mass (M > 5 ■ 10 10 Mq) galaxies; the open circles star-forming 
galaxies with e < 1; the open squares star-forming galaxies with 
1 < e < 1.3; the filled squares star-forming galaxies with e > 1.3. 
Galaxies with a mass around 10 11 Mq, corresponding in Fig.|2]to 



the characteristic magnitude M t- 



-24.7 (or L* K 



1.6 x 10 1 



Lq), have a mass-to-light ratio of \og{Mj Lk) ~ —0.2. 



of old elliptical galaxies at z > 1. In contrast, the majority 
of the star-forming galaxies have colors around 3. The color 
properties correlate with the specific SFR: The bluest ob- 
jects are also faint galaxies with the highest specific rates. 
Comparing the data from Table [1] with Fig. [4] shows that 
the counterparts are clearly bluer than the high-mass star- 
forming galaxy population, with color index lower than ~3. 
Recalling the correlation between mass and magnitude dis- 
cussed in the previous section, we note tha t the colors and 
mass -to- light ratios are also tightly related (|Bell fc de Jond 
l200lt) . There is also a strong relation, albeit with some dis- 
persion between mass and SFR, resulting in intermediate 
colors of the high-SFR galaxies. Therefore the blue and faint 
galaxies, typically characterized by high specific rates, are 
not the objects that have the highest SFRs. 

Properties like the mass-to-light ratio and color index 
are generally considered to provide information on the star 
formation history of galaxies or how they assemble their 
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Figure 5. Probability density functions of the SFR* /{SFR) (up- 
per panel) and the epoch of formation EoF, when normalized 
to the Hubble time (bottom panel), for the star- forming galaxy 
populations at z = 2.036 (dotted) and z = 0.832 (dashed). The 
counterparts of the hosts in Table [T] are superposed at arbitrary 
ordinate, but in order of decreasing redshift (from top to bottom). 



given redshift of a galaxy to its average activity since the 
object started to form its stellar populations. Both distribu- 
tions peak around unity, but the high-redshift distribution 
includes more galaxies with SFR* /{SFR) > 1, meaning 
that high-redshift galaxies are more active than present-day 
galaxies. All hosts except #2 have SFR* / (SFR) close to or 
above unity (0.85 for #9). 

The lower panel shows the epoch of formation normal- 
ized to the Hubble time. Comparing the distributions at 
the two redshifts shows that the low-redshift one is more 
extended but includes a non-negligible amount of objects 
with recent epochs of formation. All candidate hosts have 
ages within 40% of the age of the universe. We remind the 
reader that the distributions are constructed only from the 
star-forming galaxies and do not include the old, non star- 
forming galaxies. The epoch of formation indicates the epoch 
at which the galaxy was the most active and is expected to 
be quite different from the epoch of formation of the first 
stellar populations, tf orm . Table [T] shows that the epoch 
of formation of the counterparts are indeed different from 
tform- This particular time is relatively close to 1 (as de- 
fined in Table [T}, meaning that the candidate hosts include 
an early-formed stellar population, although not dominant 
since they are at the same time young objects. 

That the major fraction of the mass of most counter- 
parts was assembled in the recent times is consistently shown 
by the specific rate, the ratio SFR* /{SFR} and the epoch 
of formation. We note that the three hosts with the highest 
SFR* /{SFR) ratios (counterparts #3, 6 and 10), also have 
the highest specific rates although their SFR* and magni- 
tudes differ widely (GRB 000418, 990712 and 970508, see 
Fig. |3]) . They are among the hosts with the youngest epochs 
of formation, the bluest R—K colors, and the lowest mass-to- 
light ratios (Table [l}, pointing again to the tight correlation 
between these galaxy properties as discussed above and to 
the consistency with the GRB host galaxy properties. 



mass. They also tend to be more tightly correlated with 
the specific rate than the SFR. The SFR* /{SFR), to be 
further discussed below, is another property that strongly 
correlates with the specific rate. As already pointed out in 
Paper I, but now confirmed using the calculated observable 
properties, we note the consistency between the observed 
properties of the GRB host galaxies: If host galaxies are 
blue and faint, they are expected not to be high-SFR and 
early-formed galaxies. The properties we have found for the 
counterparts confirm this conclusion. 

3.4-2 The star-forming activity of the counterparts 

Figure [5] compares the SFR* /{SFR) and the epoch of for- 
mation of the counterparts to that of the star-forming galaxy 
population at z — 0.832 (dashed lines) and z — 2.036 (dot- 
ted lines), the highest redshift of the host sample. In each 
panel the counterparts are superposed on the distributions 
at arbitrary ordinate values, but with a decreasing redshift 
(from top to bottom). Strictly speaking only two hosts, #1 
and #6, should be compared directly to their respective dis- 
tributions. 

The top panel compares the star-forming activity at a 



3.5 Discussion 

Here, we emphasize our main result by comparing the SFR- 
to- luminosity ratios and specific SFRs of the numerical coun- 
terparts to the median values of these quantities in the cat- 
alogs they are selected from. These are displayed in Fig. HI] 
The median value of the specific SFR clearly increases with 
redshift. Note that at redshifts below z ~ 0.7, the mini- 
mal specific SFR tends to increase with decreasing redshift: 
This is due to the fact that the star formation of massive 
galaxies slows down as the redshift decreases and eventu- 
ally ceases in an increasing number of them. They therefore 
disappear from the star-forming galaxy population. Com- 
paring the specific SFRs and the SFR-to-luminosity ratio 
of the counterparts to the median values at similar redshift 
shows that the counterparts in all cases except one have 
values higher than the median. The only exception is the 
counterpart to GRB 990123 (#2). The SFR UV / {L B / L* B ) 
of this particular host is the lowest of the Chr04 sample and 
its observed spectral energy distribution was best fit by a 
star-forming Sa galaxy type, whereas all of the other hosts 
were fit by starburst templates. Th e B— band magnitu de of 
this host could however be fainter (|Bloom et al J [19991 ). giv- 
ing a higher SFR-to-luminosity ratio. Moreover the R — K 
color of this host is as blue as most of the hosts studied in 
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Figure 6. The median values (crosses) of the SFR-to-luminosity 
ratio (top) and specific SFR (bottom) for the star-forming galaxy 
populations at the redshifts of the GRB hosts. The range of val- 
ues for each catalog is shown by the length of the vertical lines. 
Diamonds refer to the numerical counterpart in each case. The 
two squares are the counterparts of GRB 000911 and 030329 hosts 
that are briefly discussed at the end of Section [3] 



iLe Floc'h etaH i|2003T ), whereas it is the reddest one of our 

sample, as seen above (see Table 1). 

It is interesting to note that ICharv et al.l (2002) esti- 
mate the specific SFRs of hosts #2, 4, 6, 7 and 8, although 
their observational study is based on extinction-corrected 
SFRs, and thus their derived specific rates are much higher 
than our results. The observed hosts #4, 6, 7, and 8 are 
among those with the highest specific SFR of the whole 
ICharv et all l|2002l ) sample. They are also found to have 
higher specific rates than local starbursts have. The numer- 
ical counterparts to #6, 7 and 8 have specific rates well 
above the median values (the first three objects on the left 
in Fig. EJ). 

The two squares in Fig. [6] show the SFR-to-luminosity 
ratio and the specific SFR of the counterparts of the hosts 
of GRB 000911 and 030329, discussed at the beginning of 
this section. These counterparts were searched for in cat- 
alogs at slightly different redshifts from the host redshifts, 
z — 1 and z — 0.133, respectively. The counterpart of GRB 
000911 (AX = 0.3% and AY = 5.2%) is a low-mass, young 
object (M = 1.1 xlO 9 M , l-EoF/t H(z) =0.37) with a high 
specific rate (e=1.87), blue color R — K = 2.12 and a mass- 
to- light ratio M/Lk = 0.36. The closest counterpart that 



we find for the host of GRB 030329 has a lower SFR-to- 
luminosity ratio than the observed estimate (AX = 1.7% 
but AT" ~ 19%). The candidate host is a low-mass ob- 
ject (M = 1.4 x 10 9 M ), with blue color R - K = 2.2, 
a specific rate of e = 0.7 and an epoch of formation of 
1 — EoF/tjj( z )— 0-45. Compared to the overall population 
these two candidate hosts have SFR-to-luminosity ratios and 
specific SFRs well above or similar to the median values at 
the same redshifts. 

The fact that the Chr04 sample is magnitude limited 
and includes bright hosts may explain why the numerical 
counterparts are not found among the objects with the high- 
est possible specific SFRs. The counterpart of GRB 970508 
(#6) is an example. From HST observations and other stud- 
ies, thi s host may be descr ibed as a blue compact dwarf 
galaxy l|Fruchter et al.|[2000t ). Could this host be a prototype 
of the general host GRB galaxy population? Interestingly, 
the counterpart of this host could be found in the sample 
of com pact galaxies at z ~ 0.7 analyzed in iGuzman et al.l 
(1997). The counterpart #6 (z = 0.83) with an absolute 
magnitude Mb = —18.2, a mass close to 10 9 Mq and a spe- 
cific rate of e ~ 1.8 falls into the fainte st and highest specific 



SFR population of compact galaxies in Guzman et al 



(see their figures 5 and 7). Moreover. ISollerman et al 



1997) 



2005) 



discuss three GRB host galaxies and show that they have 
similar properties as a sample of compact blue galaxies in 
the local universe. 



4 CONCLUSIONS 

In this paper, we have discussed cosmological galaxy prop- 
erties and the properties of host galaxies of GRBs in partic- 
ular, using fully hydrodynamic simulations of galaxy for- 
mation and the stellar population synthesis (SPS) code, 
Starburst99 , to infer observable properties. An important 
feature of the numerical procedure is that the star forma- 
tion history of galaxies entering the SPS code, comes di- 
rectly from the simulation, with each stellar population con- 
tained in a given galaxy treated as a homogeneous popula- 
tion, that forms instantaneously. We identify objects in the 
simulation that have optical star formation rate, SFRuv, 
and SFR-to-luminosity ratio, SFRuv /(Lb /L%), similar to 
those estimate d in a well-defined sam ple of ten observed 
host galaxies (jChristensen et all |2004| ) , the only available 
homogeneous sample focusing on this ratio. Each numerical 
counterpart is selected from a simulated catalog at the same 
redshift as the corresponding observed host, and is charac- 
terized by a number of properties: The SFR* , mass, specific 
SFR, SFR* /(SFR), epoch of formation of the first stellar 
population, and epoch of formation of the galaxy object. 
These are obtained directly from the simulation. In addition, 
B and K— band luminosities, i?-band apparent magnitude, 
R — K color, and mass-to-light ratios are obtained from the 
SPS or combination of both results. It should be emphasized 
that some of these properties (e.g. mass and SFR*) are esti- 
mated directly from the simulations, making their definition 
inherently different from those commonly adopted in obser- 
vations. 

The IChristensen et all l)2004h sample includes host 
galaxies with redshifts in the range 0.43 < z < 2.03. The 
sample is magnitude-limited (R < 25.3), with estimated ab- 
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solute B— band magnitudes between —21.4 and —18.1. Our 
counterpart hosts are low-mass galaxies (M < 4 • 10 10 M©), 
with low mass-to-light ratios (M/Lb around 0.5); most of 
them are blue (R—K < 2.9) and young galaxies, with epochs 
of formation (or ages), within 40% of the age of the universe 
at the different redshifts. Although the SFR* of the coun- 
terparts varies between ~ 0.4 and 8 M© yr -1 , the specific 
SFR is equal to or higher than the median values estimated 
for the different catalogs, with the lowest value being e ~ 1. 
Because of its strong correlation with the specific rate, the 
SFR* /{SFR) also has high values, around unity or higher. 
To outline the consistency of such an ensemble of properties, 
we compare the counterparts to the overall galaxy popula- 
tion at intermediate redshift and discuss the strong rela- 
tionships between the specific SFR and quantities known to 
reflect the star formation history of galaxies, i.e. color and 
mass-to-light ratio. Indeed, the bluest objects are also faint 
galaxies with the highest specific rates. 

Our identification of simulated galaxies with observed 
hosts has some limitations, both from observational and nu- 
merical points of view. Some of the observed hosts may, in 
other observational studies, be found to have slightly differ- 
ent SFRs or magnitudes, making their SFR-to-luminosity 
ratios uncertain. Although the general agreement between 
the simulation and the observations is fairly good, it should 
be kept in mind that the moderate resolution and some- 
what limited number of physical processes included in the 
simulation may affect the relative number of galaxies in 
each sub-population. For instance, these simulations only 
account for a limited number of "extreme" -type galaxies, 
such as extreme starbursts or massive star-forming galax- 
ies. F urthermore, we do not in clude the effects of dust and, 
as in IChristensen et all (|2004f ) , we only consider extincted 
SFRuv in the estimate of the SFRuv / (Lb/L* b ) ratio. That 
may not be a serious drawback since our focus is on low-mass 
galaxies that on average suffer less attenuation than the mas- 
sive ones. In addition, the amount of dust h as been shown to 
be lim ited in most host galaxies of GRBs (|Le Floc'h et al.l 
2006). The main limitation of this study is the wide red- 
shift range of the 10 observed host galaxies, with only a few 
galaxies in each narrow redshift bin. The inferences made 
therefrom need to be confirmed once the sample size has 
grown by a factor of 5-10, with a number of hosts at sim- 
ilar redshifts. We should expect to see still higher specific 
rates in low-mass objects, as the simulated catalogs tend to 
show. The evolutionary effects displayed in Fig. [6] support 
this conclusion. 

Comparing an expanded host galaxy sample based on 
Swift data with the simulated catalogs and other well-known 
galaxy sub-populations, such as compact blue galaxies may 
be extremely useful. Hosts galaxies should provide a clearer 
view into the formation and evolution of galaxies and the 
role of different sub-populations therein. In particular, they 
will help in investigating the faint end of the galaxy lumi- 
nosity function. Our results, obtained using a numerical ap- 
proach, are consistent with and confirm the picture of GRB 
host galaxies that has emerged lately: They tend to have low 
mass, be blue in color and have relatively high specific star 
formation rates. Host galaxies may then belong to the high 
specific SFR galaxy population, rather than the high SFR 
population. High-resolution simulations are required to de- 
termine, in a more quantitative way, the contribution of the 



host galaxies to the overall population, and whether hosts 
are a part of the average normal star-forming galaxy pop- 
ulation or a sub-population of this one, with blue and very 
low-luminosity objects. 
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